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O I The functional relations of the transfer matrices of fusion hierachies for six- 

, and eight-vertex models with open boundary conditions have been presented in 

this paper. We have shown the su(2) fusion rule for the models with more general 
reflection boundary conditions, which are represented by off-diagonal reflection 
^ . matrices. Also we have discussed some physics properties which are determined by 

. the functional relations. Finally the intertwining relation between the reflection K 

O ■ matrices for the vertex and SOS models is discussed. 
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Oh! 

^ I Great progress has been made on the study of integrable models in statistical me- 
chanics and quantum field theory. It is very clear that the important methematical 
^ : structure that ensures the exactly solvablity of these models is governed by the Yang- 
. Baxter equations if the models sit on square lattices with periodic boundary conditions. 
Recently there has been interest in studying integrable systems with open boundary con- 
ditions. The open boundary conditions are described by boundary reflection matrices 
satisfying reflection equations [Q, ^ (boundary Yang-Baxter equations), which ensures 
the exactly solvablity of the models with the open boundary conditions together with the 
Yang-Baxter equations. The boundary reflection matrices have been found for many in- 
tegrable systems, in particular for the six- vertex and eight-vertex models in [^, ^, ^ (also 
see 1^ 1^ for related works). The eigenvalues of the transfer matrices have been solved 
for the six- vertex model or Ai'^ invariant chain 0, |^, |13| (also see |]16| for related work), 
Izergin-Korepin vertex model or A'^^ invariant chain |lO], [/g(sp/(2, l))-invariant t-J 
model|Tl|, [1^, A^^) invariant chain ^ and Aj^^ invariant chain [Q. These exact solutions 
are constructed by the Bethe ansatz only for diagonal reflection matrices. The integrable 
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systems with non-diagonal reflection matrices are more difficult to solve and not many 
results have been obtained for such systems. 



The fusion procedure has been shown very useful in studying two-dimensional inte- 
grable models. Similar to the method to fuse the R matrix of the Yang-Baxter equation 
flSlI , the fusion procedure for the reflection matrix K of the reflection equation has been 
presented in |]T5|]. The fused R- and i^'-matrices generate some new integrable models 



with the open boundaries based on the elementary model. The corresponding fused 
transfer matrices of fusion hierarchies are related through the functional relations, which 
can be shown by the fusion. In [El] the functional relations of transfer matrices for the 



six-vertex model with diagonal reflection matrices have been shown and thus the finite 
size corrections to the transfer matrices have been obtained by solving the functional 
relations. The functional relations for the six-vertex model and eight-vertex model with 
non-diagonal reflection matrices have not been found even though the fusion procedure 
was described generally before. 

In this paper we study the six-vertex model and eight-vertex model with the non- 
diagonal reflection matrices found in 0, ^ . We find the functional relations among the 
eigenvalues of the transfer matrices for whole fusion hierarchies of the models. This sheds 
light on the solution of the models with boundaries via the most general non-diagonal 
reflection matrices. In order to show the derivation of the functional relations of the 
six-vertex model and eight- vertex model with non-diagonal reflection matrices clearly we 
carry out the fusion of the models. Then the intertwining relation between vertex and 
SOS models for the boundary reflection K matrices is also discussed and the reflection 
equations for the SOS models have been explained. From the functional relations we 
derive the crossing unitary conditions of bulk and surface free energies. For the six- 
vertex model the crossing unitary condition has been argued in terms of quantum field 
theory in PU[ and however it is not very clear for the {of f -critical) eight-vertex model. 



Our crossing unitary conditions of bulk and surface free energies must be consistent with 
the argument given in for the six- vertex model. We believe the crossing unitary 



conditions for the off-critical model is a new result. Specially, some surface critical 
behaviour can be studied by solving the crossing unitary conditions. Also we present the 
eigen-spectra of the transfer matrix of the eight- vertex model with the diagonal reflection 
matrices, which solves the functional relations. The eigen-spectra has not been found 
before. 

In section 1 we describe the fusion of the R and K matrices obviously. Then the 
functional equations are presented in section 2. For clarity all proofs are given in Ap- 
pendices. In section 3 we present the functional equations for the six-vertex model and 
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eight-vertex model with some known boundary reflection matrices. In section 4 the in- 
tertwining relation between the vertex boundary K matrix and the SOS face boundary 
K matrix is discussed. Finally we conclude with a brief discussion at the crossing unitary 
conditions for the surface free energies and the eigen-spectra of the transfer matrix and 
discuss some physics consequence. 

1. Some expressions for fusion 

Let us begin by reviewing the fusion procedure for a given R matrix obeying the Yang- 
Baxter equation 

R^\u)R^\u + v)R'\v) = R^\v)R^\u + v)R^\u) (1.1) 

where m, v are the spectral parameters. Suppose that i?-matrix acting on C^x satisfies 
the following properties, the PT symmetry 

Pr2B'\u)P^2 =Piu)= R^\u) (1.2) 
where P12 is the permutation matrix, which can be represented by 

i?i2(o) = W)Y"p^2 , (1.3) 

the unitary condition and the crossing unitary condition 

R^\u)R^\-u) = p{u) (1.4) 

R^^ (u) R^^ {-u - 2A) = p{u) (1.5) 

where ti denotes the transposition in the first space and A is crossing parameter. These 
p{u),p{u) are some scalar w-dependent functions. 

Fusion is the idea to build up new solutions of the Yang-Baxter equation. Define the 
projector 

Y+ = i(pi'f + ... + pP-i'P + /)... (pi'2 ^ /) (1.6) 
^ pi 

the new solutions can be given by [|18 



R^iu) = i;+%,g)(M-gA + A)---i?(p,2)(M- A)%,i)(M)y+ (1.7) 
R^p^j){u) = Y+R'^\u)R^'^{u + \)---RP'^{u + p\-\)Yp+ (1.8) 

where R(p^q){u) acts on C^+^ x C^"'^. These fused R matrices satisfy the fused unitary 
condition and the fused crossing unitary condition 

= pMypy, (1-9) 



where 

~p,» = iiii~pU+ik~j)x) . (1.11) 

Note that the superscripts 1, 2 mean R}lq)iu) e CP+^ ® C'?+^ We suppress them if there 
is no confusion. 

The reflection matrices K-{u), Kj^{u) satisfies the reflection equations [|T], Q 

R^\u - v)K\{u)R^\u + v)Kl{v) 

= Kl{v)R^\u + v)K\{u)R^\u - v) (1.12) 

R^^{y-u) Kl{u)R'\^-u-2\) Kl{v) 

=Kl {v)R^\^-u-2\) Kl {u)R^\v-u) . (1.13) 

Similar to the fusion of the R matrix the reflection matrices K_{u), K^{u) acting on 
space can be fused to give new solutions of the reflection equations. Following |15| 
the new solutions of the reflection matrices can be expressed by 

K^J'\u) = p{u\q) Y, [Kl{u)][R'^''-\2u + \)K'L~\u + A)] 

[R'i^'i-^{2u + 2A)i?^-i'^-=^(2M + 'i\)Kl-^{u + 2A)] ■ • • 
[i?«'i(2M + qX- X)R''-^'\2u + gA) ■ ■ ■ 

R^'\2u + 2gA - 3X)K^_{u + qX - A)] Yq (1.14) 

K^_^\u) = Yq [Kl{u) k-^'\~2u - 3A) ■ • • k'\-2u - qX) k'''{-2u - gA - A)] ■ ■ ■ 
[Kl{u + qX- 3A) k'^{-2u - 2qX + 3A) k'^{-2u - 2qX + 2A)] 
\Kl{u + qX- 2A) k'\-2u - 2qX + X)][KI{u + qX - X)]Yg (1.15) 

where K^^\u) , K^^\u) act on O+i, k'^ Jr''^ and 

p(u\q + l) = p{u\q)/pq^i{2u + qX) 

p{u\l) = 1, g = l,2,--- . (1.16) 

The fused reflection equations become 

R(P,g)iu - v)K^l\u)R(q^p){u + v- X+pX)K^J'\v) 

= K^^\v)R^q^p){u + v-X+pX)K^:'\u)R^p^q){u-v) (1.17) 

tl t2 

i?(,,p)(i; - u) irf (u) %,,)(^ - - A - pA) Ki'^ (v) 

= Ki'^ {v)R(p,q){^ -u-X-pX) K^^ {u)Ri^q,p){v - u) . (1.18) 

4 



There is an automorphism between the relation for K^^^ and i^!!'"' [|, 



K^:^\u)/p{u\p-l) = K^l\-u-p\) 
K^^\u) / p{u\q - I) = K^l\-u-qX) 

and the automorphism for i?-matrices 

^(g,p)(^) = + gA - p\) 

which can be seen directly from the fusion procedure. 



(1.19) 
(1.20) 

(1.21) 



The fused Yang-Baxter equation and the fused reflection equation for R and K ma- 
trices guarantee the following commuting families 



u 



J'(P:6)( 







where 



T(^''«)(n) = tr( ki'^\u)Ui^,,,){u)K^'\u)U^,,^){u + qX - \) 



(1.22) 



;i.23) 



U 



(P,.)1'«J = + - \mtp)^u + IpA - iA) ■ ■ ■ , Rll^{u + ipA - \\) (1.24) 

-,Af,c 



C>(P,.)(^) = - IPA + iA) ■ ■■Rl],){u - ipA + \\) R\;^^){u - \p\ + \\) (1.25) 

and & = 1,2,---. The proof of (|1.22|) can be done similarly to the study of the 
unfused six-vertex model given by Sklyanin and is briefly described in Appendix A. 
using graph representation. 

2. Functional equations for general reflection matrices 

The last section has shown that fused models can be built up by fusion of the elementary R 
and K matrices. The fused models are integrable systems because we have the commuting 
families ( |1.22|) . In fact, the these fused transfer matrices are related each other by the 
groups of functional relations, which can be constructed by fusion. For the periodic 
boundary condition it was already well known that the transfer matrix 



u) 



tr C/(p,g)(n) 



commutes 



and the functional relations 



I^(p,g)(^) , T(P''?')(y) 



u 



(2.1) 



(2.2) 



(2.3) 



^some arbitrary parameters in K)^ are interchanged, e.g. (^+,/i+,j^+) ^ for the six 

vertex model (3.3) 
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can be constructed |T^, where T^^''^\u) = / G C'^'^^ and the w-dependent function 
is generated from the antisymmetric fusion of the modeL These relations, in fact, are the 
su(2) fusion rule and show the implication of the eigenvalues of transfer matrices of all 
fusion hierarchies. They have been applied successfully to derive the eigenvalues and the 



relevant Bethe ansatz equations of the fused or the elementary transfer matrices [T^, ^0 
Also they could be solved directly in the thermodynamic limit to yield the finite size 
corrections to the eigenvalues of the transfer matrices ETI, E^l, EH 



For the models with open reflection boundary condition it has been noticed that the 
functional relations can be applied to obtain the eigenvalues of fused transfer matrices 



T3| , |T5[ and to find the finite size corrections to the eigenvalues of transfer matrices ^4 
These studies have concerned only models with diagonal reflections. For the off-diagonal 
reflection boundary solvable models we still expect the su(2) fusion rule. This is shown 
in this section. 

Let us first prepare some notations as follows. 

N 



i=l 
N 



I+{u\q) := K,~i?;i7_\)(gA - 4A - 2u)Rlf-}^^{q\ - 3A - 2^)1^2- e (2.4) 

"^""'^^ ■ p,,i(2a-gA+2X)p^i,i(2a-gX+A) 

where Y-f = |(1 — P^'^) is an antisymmetric projector. The boundary dependent nota- 
tions are defined by 

iu+{u) := kl{u + A)P^'^(-2m - 3A) kl{u)Y.^ 

uj_{u) := Y^Kl{u)R^''^{2u + X)K'i{u + X)Y^ (2.5) 

The product of all of these notations is expressed by 

fiu) = uj_{u)uj+{u) {^ci)l{u)ctf_{u)j ® {l+{u\q)I^{u\q)^ E ^^^.C^^^ ® , (2.6) 

which is clearly a g-independent function if we have the relation 

/+(«|g)J_(n|g) = Jpi;i(2n + A) , (2.7) 

where the identity matrix depends on the fusion level g or / G C^. The above relation 
is just right for the six- and eight-vertex models, which will be seen in the next 
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section. Moreover (p and / are the diagonal matrices and so is f^{u), which is written in 
the form 



r{u) = i-nu) 

nu) = uj.{u)uj+{u)cl>l{u)cl>Uu)lMq)I-{n\q) 



(2.8) 



and I G ®^iC^^^ ® C^. Now we can express the functional relations for the six- and 
eight- vertex models with open reflection boundary conditions in the following theorem. 



Theorem 1 (sn(2) Fusion Hierarchy) Let us define 

F, = F{u + qX) , T(P'°) = I e ®l,C^+' 
with the transfer matrices 

given by ( {1.23^ ). Then the su(2) fusion hierarchy follows as 



rp{q) T-i(l) _ I fP rjiil 

-'■0 -'-a — -^0 "T Jn~l -L 



(9-1) 



q = l,2, 



(2.9) 



if the reflection K_,K+ matrices satisfy the reflection equations (\1.12^ )- ^.l^ ) 



The theorem is proved in Appendix B. Then following the standard procedure given by 
21| , |25[| by introducing the inversion identity hierarchy (or y-system) 



JO) 







(2.10) 
(2.11) 



fc=0 



the functional equations for t can be derived, which are called thermodynamic Bethe 



ansatz (TBA) equations in |21|, Q . To see this let us consider the triple transfer matrices 

jiCp.?)^ rjp{p,q-l) J^(p,l]^ _ ^ rjp{p,q) j^{p,l)^ j,{p,q-l) 



Inserting the functional relations (|2.9|) into the terms in parentheses this equation gives 
new functional equations, 



-^1 



k=0 



(2.12) 



This leads to the following theorem. 
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Theorem 2 (sn(2) TBA) The inversion identity hierarchy satisfies the following 
su{2) TBA equations 

ti'h['^ = (/ + + t['-'^) (2.13) 

for any reflection K_, matrices satisfying the reflection equations (j7T7^-([7TT^, where 



3. Examples 

We have shown that the functional relations hold for the six- and eight-vertex models 
with any K reflection matrices satisfying the reflection equations ( p..l2| )-( prT5 ). The su{2) 
fusion rule depends on the boundaries though the function /^(m) defined in (p.8| ). In this 
section we calculate explicitly the function for the six- and eight-vertex models. Here we 
use the general K matrices presented in [Q and 

3.1 Six vertex model 

The six vertex model is described by the following R matrix, 



R{u) = 
a{u) = sm{u + A) 



fa{u) 
b{u) c{u) 
c(n) b{u) 

V 

b{u) = sin(M) 



\ 



a{u)J 
c(u) 



(3.1) 



sin(A) 



and the boundary reflection matrices are 
Kju) = 



sin(^4. — u — X) z/_(_ sin(2u + 2A) 
/i+ sin(2M + 2A) sm(^+ + u + A) 
sin(^_ + u) /i_ sm(2u) A 
sin(2u) sin(^„ ~ ^) / 



(3.2) 
(3.3) 



where ^±,/i±, i'± are arbitrary parameters. 

The unitary conditions can be checked by inserting the R matrix 
(11. 51) and we obtain 



Ah into (O) and 



p{u) = sin^(u) — sin^(A) ; p{u) = sin^(A) — sin^(u + A) = —h{u)a{u + A) (3.4) 
The quantities in ( |2.4| ) are given by calculation to be 



cc;+(m) = — sin(2u+4A) ^sin(,^++M + A) sin(,^+— m — A) — /i+z/+ sin^(2n + 2A) 
uj^{u) = sin(2M) ^sin(^„+M + A) sin(^_ — n — A) — /i_z/_ sin^(2u + 2A)^ (3.5) 
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p-1 



0+H= n[-pi,i(«-fcA+ipA-iA)]^ 

k=0 

= n [-pi,i(« +kx- bx + 

k=0 
q-1 

k=l 

i~{u\q) = - f[[-PiA2u + 2X-kX)]-K 

k=l 



(3.6) 



(3.7) 



It is obvious that the relation (|2.7|) is satisfied. 



The fusion procedure for the R matrix brings the extra zeros to the fused R matrices. 
The number of the zeros are only dependent on the fusion level p, q. Some zeros of 
R(^p^q){u) can be removed from U(p^q){u) and U (^p^q){u + qX) in ( |1.24| ) and ( |1.25| ) by the 
replacements 



{9:P) 



.U 



Ri,A^)/'W]\Hu + 3X-kX)] 

j=0 k=0 

RM{u)/lili[b{u + jX-kX)] 
fc=oi=i 



(3. 



(3.9) 



Rip,q){u) 

The replacement only changes 4'±{u) to a simpler form, 

(Pl(u) = (P'Liu) = [b{u - yX + iA)a(n + ipA + |A)]^ 
Inserting all of the quantities given in (|3.5|) - (|3.9|) into the definition ( p.8|) we have 

r{u) = uj+{u)uj^{u) ^ I I 3.10 

b{2u + X)a{2u + 2A) 

For the simple case of diagonal K matrices the fusion rule ( p.9|) with (|3.10|) coincides 
with that obtained in 123 . 



3.2 Eight-vertex model 

The eight- vertex model is the generalization of the six- vertex model The R matrix 
for the eight-vertex model is given by 

fa{u) d{u)\ 

b{u) c(n) 

c{u) b{u) 

\d{u) a{u) / 



R{u) 



(3.11) 



These non-zero elements are given by 

a{u) = H{u + A)e(M)e(A) , b{u) = e{u + X)H{u)e{X) , 
c{u) = e{u + X)e{u)H{X) , d{u) = H{u + X)H(u)H{X) 



(3.12) 



9 



where H{u), Q{u) are the theta functions. The relevant boundary reflection matrices are 
given by g] 

/ H{C-+n)e{^.-u) /._A:V2sn(2^z)sn2(n)//(n;C_)\ 
-^''> [-f,^k-y'sn{2u)/f{u;C-) Q{^. + u)H{C. - u) ) ^^'^^ 

K^{u) = (3.14) 
/ H{^+-u- X)Q{^+ + u + X) ^+A:-V2sn(2n + 2A)//(n + A;e+)\ 

\-^+A;i/2sn(2u + 2A)sn2(n + A)//(n + A;e+) e{C+ - u - X)H{^+ + u + X) ) 

with 

/(n; = [1 - sn^H sn2(0]/[e(e + «)e(e - u)] , 

and 

sn(u) = k-^/^H{u)/Q{u) 

where /i± are arbitrary parameters. The scalar functions in the unitary relations ( |1.4D 
and ( |1.5| ) for the eight vertex models are 

p[u) = h{\ + u)h{\-u) (3.15) 
p{u) = h{u)h{u + 2X) (3.16) 
h{u) = H{u)e{u)e{0) (3.17) 

The expressions of the quantities in ( ^^ with ( p.l6[ ) are still correct except the uj±{u), 
or 

uj_{u) = h{^_ +U + A)/i(^_ - M - X)H{2u)e{2u + 2A)/e(0) (3.18) 

+K4u)2iK4u + X)2iH{2u + 2X)H\u + X)e{2u)e{0)/e\u) 
u;^{u) = -h{^+ -u- X)h{^+ + U + X)H{2u + 4A)e(2u + 2A)/e(0) (3.19) 

-K+{u + X)i2K+{u)i2H{2u + 2X)H\u + A)0(2m + 4A)e(O)/0=^(M + 2A) 

p-i 



u 



^[-~p,,{u-kX + \pX~\\)]' 

p-1 



[u 



^[-p,^,{u + k\-yx + \\)f (3.20) 

fc=0 
9-1 

I+{u\q)= n[-Pi,i(2w + A-/cA)] 

k=l 

q 

I.{u\q) = - n[-Pi,i(2M + 2A- A;A)]-^ (3.21) 
fe=i 



From (|3.21| ) it can be seen that (|2^ ) is satisfied. 



Similar to the six-vertex model discussed in the last subsection, removing the zeros 
generated by fusion from U(j)^q^{u) and t/(p^g)(M + qX) in (|1.24] ) and ( |1.25| ) 



g-lp-2 

R(,,p){u) Ri,,p){u)/ n ElHu + jX - kX)] 

j=0 k=0 
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7-1 p-1 

fc=oi=i 

from the transfer matrix ( |1.23| ), the quantities (p±{u) become 



(3.22) 



(Pliu) = [h{u - yx + ix)h{u + yx + + \)Y 

0P(m) = [h{u-yX + lX)h{u + yX + lX + X)f . (3.23) 
Inserting all of the quantities given in ( p.5| )-( |3l9D into the definition ( p.8| ) we have 



[h{u - yx + ^X)h{u + yx + |A + A)]^^ 
h{2u + X)h{2u + 3A) 



(3.24) 



4. SOS models and intertwiner 

Intertwining relations are the correspondence between the i?-matrix (the Boltzmann 
weights of the vertex model) and VT-matrix (the Boltzmann weights of the relevant face 
model) |29| . Baxter introduced the intertwining relation for the eight- vertex and the 
relevant face model to solve the eigenvalues of transfer matrices of the eight- vertex model. 
The related face model is called the solid-on-solid (SOS) model. The similar formulation 
can be generalized to the boundary K matrices (boundary Boltzmann weights). This 
however has not been done before. We show briefly in this section the correspondence 
between the boundary K matrices of the eight-vertex and the relevant SOS models. For 
clarity the corresponding graph representations are briefly given in appendix C. 

At first let us review Baxter's intertwining relation between the eight-vertex and the 
SOS model. Choose arbitrary constants and integers a, 6 G Z and set the nonzero 
intertwiners by 

^aAu)={lY^"~'''\]/4s^) ife = 6-a = ±l (4.1) 
\^ B(s + a — eu) J * 

where the factors S{b) is given by 

S{b) = h{wo + bX) (4.2) 



The R matrix is given by ( |3.11| ) with (|3.12| ). Then Baxter's intertwining relation is given 
by 



R{U - V) ipdM ® Vc,b{^) = 'PdA'") ® fafiiu) W\ 



fd 


c 










u - 


..) 


\ a 


h 







(4.3) 



where W 
relation 



fd 


c 










\ a 


b 





u are the face weights of the SOS model satisfying the star-triangle 



7 9 
a b 



u\W 



e d 
,/ 9 

11 



V ]W\ 



d c 
9 b 



v — u 



e g 



v-u \ W\ 



9 c 
a b 



9 \f « 

for any a, b, c, d, e, f. The nonzero weights are given by 



V \W 



' e d 
.9 c 



u 



(4.4) 



a a ± 1 

a±l a±2 

a a ± 1 
a ^ 1 a 

a a ± 1 
a ± 1 a 



u 



u 



h{X + u) 



S{a + l)S{a-l) h{u) 



u = 



h{wo + aA =F m) 
h{wQ + aX) 



h{X) 



(4.5) 



where a e Z and wq = ^{s~^ + s )X — K. These face weights satisfy the following unitary 
condition 



u 


c 




a' 










\ a 


b 




b 



c 

and the crossing unitary condition 

c 

where the scalar function is given by 



/ a 


d 




d 








-u-x)w[l 


a' 


u 


\b 


c 







-U = p{u)5a 



\ S{a)S{c) 



(4.6) 



(4.7) 



p{u) 



h{X - u) h{X + u) 



h{X) h{X) ' 
The intertwining relation for the boundaries is set out by 

K_{u)ipb^a{u) = ^ipb^c{-u)K-(b 



(4.^ 



u 



t S{c) ( c 



u 



(4.9) 



where K_ ( b 



u 1 and Kj^ { b 



uj are the reflection matrices {K matrices) of the SOS 

model. It is obvious that the face K matrix elements are nonzero only for |6 — a| = 1 
and \b — c\ — 1. The nonzero elements are found to be 



b 



= ^-l(u + X)K+{u)^b,c{-^-X)^ (4.10) 
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The intertwiners are dependent on two arbitrary parameters and thus (pa,a+i{u) 
and (pa,a-i{u) can be treated as independent vectors. We therefore can define the inverse 
vectors of them by 



J2 'Pa,b{'^)i'Paiiu)j = 



b=a±l 



1=1,2 



where 

Valiu) = S{b)(^ee{s-' + a + eu), -eH{s-' + a + eu)^ /deta['^{u)] 
with e = 6 — a = ±l and 

2H{x - K)e{x - K)H{y)e{y) 



{4.11] 



(4.12) 



with X = h(s' 



deta[ip{u)] -- 
s+) + aX and y 



H{K)Q{K) 
s+) +u. 



Applying the intertwining relations ( [4.3| ) and ( ^.9| ) to the refiection equations for the 
vertex model we are able to obtain the refiection equations for the SOS model, 



f,9 



a b 
.9 c 



-UK- 

f,9 

'a g 
b c 



u — V j ( g 



f 

v — u \K 



c 
f 

a b 

a f 



M + ( g 



f 



u]W\ 



f 

e 

a g 



(4.13) 



f,9 



f 



v]W 



a b 
9 f 



u]W\ 



a g 
d f 

-2X-U-V ] K+ 



f 



d e 
d 



u — V 



9 



u]W\ 



a g 
d e 



J S{d)S{g) 

S{f)S{a) 



v — u 



S{b)S{g) ■ 



This SOS analogue of the refiection equations, directly following from the refiection 
equations of the vertex model, has not been given before Similar to the vertex 

model, we have the commuting transfer matrix V{u) defined by the following elements 



{a\V{u)\b) 



{co,--,cn} 

xWi 



fao 



Co 



u 



\bo 



n 

fc=0 



W\ 



Ck Cfc+1 



Jk+1 



U + Vk 



U-Vk 



KJcm 



U 



which satisfies 



[V{u), V{v)]=Q 



(4.15) 



(4.16) 



where a = {qq, ai, • ■ ■ , a^}, b = {bo, &i, ■ ■ ■ , &7v} and these Vk are some arbitrary parame- 
ters. The proof of ( [4.16|) can be done similarly to the vertex models. Also the su(2) type 
fusion rule in theorem (|T]) still works for the SOS models with open refiection boundaries. 
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Theorem 3 (sn(2) Fusion Hierarchy) Let us define that 



are the fused transfer matrices of the SOS model with the fusion levels p and q. Then the 
su(2) fusion hierarchy follows that 



yi^) yil) ^ + fP yi^ 



(9-1) 



(4.17) 



if the reflection K^, matrices satisfy the reflection equations l \4-13i) -( \4-i4 )- The func- 
tion f^ = f^{u + qX) is given by the antisymmetric fusion. 



The proof can be done similarly to the case of vertex models. In appendix C we 
give the proof for the simpler case of g = 1 and p = 1. It is interesting to notice that 
the functional relations ( [4.17| ) of the SOS model with the open boundaries and periodic 
boundaries have the same form except the w-dependent function fg_i being different. 

5. Discussion 

In the previous sections we have studied the fusion hierarchies of the six and eight vertex 
models with the open boundaries. The functional relations of the fused transfer matrices 
of the models with the non-diagonal reflection matrices have been obtained. It has been 
shown that the functional relations can be solved to obtained the finite size corrections 
to transfer matrices for the six- vertex model with the diagonal refiection matrices ||2^. It 



should be interesting in the near future to solve the fusion hierarchy and thermodynamic 
Bethe ansatz equations for the six-vertex model with the general non-diagonal refiection 
matrices by similar techniques. In this way it should be possible to obtain the central 
charges and scaling dimensions in terms of Rogers dilogarithms and their analytic con- 
tinuations and hence completely elucidate the critical behaviour of the model and their 
fusion hierarchies. 

The free energies and the eigen-spectra of the transfer matrices of the models with 
the open boundaries can be determined by the functional relations. To show this the 
discussion can be divided into a number of points. 

5.1. Crossing unitary condition for boundaries: The crossing symmetry for 
boundary K matrices have been expressed in Using the crossing symmetries and 
the unitary conditions of the boundary K matrices the surface free energy fs can be 
fixed. In the argument is given in terms of quantum field theory. According to this 
argument it is not very clear how to express the crossing and unitary relation of the 
surface free energy for an off-critical model. However, like the R matrix the crossing 
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symmetry of the K matrix can be given by the fusion procedure in statistical mechanics. 
To see this let us first consider the crossing symmetry of the R matrix. 



The antisymmetric fusion of the R matrix 

YiR''\u - X)R''\u)Y,- = -p,^,{u - A) = p,^,{u) (5.1) 

is equivalent to the crossing unitary condition ( |1.5| ). This relation leads to the crossing 
unitary condition for the bulk free energy 

/,(w)/,(n + A) = [pi,iHr (5.2) 

Similarly, the crossing unitary condition of the K matrices is related to the antisymmetric 
fusion of the K matrices (see ( |2.5| )). More precisely the factors u± determine the crossing 
unitary conditions of the K± matrices. Thus the surface free energies are determined 
by 

f^^-)f^^- + >^)= h(2u-X)h[2u + X) 

The denominator h{2u — X)h{2u + A) or sin(2M — A) sin(2M + A) in the function f^{u — A) 
can be removed by normalization of the transfer matrix T^^'^\ Thus the functions f^{u) 
in ( |3.1U| ) or ( |3.24| ) contain the unitary conditions and crossing symmetries of the bulk R 
and the boundary K matrices. 



The bulk free energy of the eight- vertex model has been obtained by Baxter . By 
the similar method we can solve the the ( |5.3|) to obtain the surface free energy for the off- 
critical model. Therefore we should have the surface specific heats and the related critical 
exponents. This study does not rely on the exact diagonalization of the transfer matrix 
of the model with the non-diagonal refiection matrices, which has not yet solved. Thus 
we have presented a possible method to find the crossing unitary conditions for solving 
the surface free energy of an exact solvable model with open boundary conditions. 

5.2. Ansatz of eigen-spectra: We have found the functional relations of the transfer 
matrices for all fusion hierarchies of the six- and eight- vertex models with the most general 
reflection matrices. The su(2) fusion rule still hold for the models with any K matrices 
satisfying the reflection equations. It is not very clear to extract the exact eigen-spectra 
of transfer matrices from the functional relations for the models with general reflection 
matrices. However we may solve the eigen-spectra of transfer matrix from the functional 
relations for the eight- vertex model with the diagonal reflection matrices, which, in fact, 
has not been given before . 

For the special case of diagonal K matrices, the transfer matrix T(u) = T^^'^\u) 
of the six vertex model has been diagonalized using algebraic Bethe ansatz and the 
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eigen-spectra is given by 



T{u) = (Ai(n) + A2(M))/(7(n) (5.4) 

Ai(u) = uJi{u)(p{u + \)Q{u- X)/Q{u) (5.5) 

A2(m) = uj2{u)(l){u)Q{u + \)/Q{u) (5.6) 

g{u) = sin(2M + A) (5.7) 



where 



M 

Qi.'^) = Yl {sm{u - Vm) sm{u + Vm + X)} (5.8) 

m=l 

0(u) = sin^^iu) (5.9) 

uj2{u) = ti;i(— M — A) = sin(,^+ + M + A) sin(^_ — M — A) sin(2M) . (5.10) 

These vi,V2, ■ ■ ■ ,vm satisfy the Bethe ansatz equations 

T{vj) = . (5.11) 

It is easy to see that 

A,{u + X)A2{u)=f\u). (5.12) 



Thus ansatz ( ^.4|) solves the functional relations (|2.9| ). If we use the same ansatz ( ^.4|) 
with all trig-functions sin replaced with elliptic function h, 

M 

Q{'^) = Yi {h{u - Vm)h{u + Vm + X)} (5.13) 

m=l 

cp(u) = h'^'^iu) (5.14) 
U!2{u) = uJi{—u — \) = h{^^+u + \)h{C,_—u — X) 

xH{2u)e{2u + 2\) , (5.15) 

g{u) = h{2u + \)e{0) (5.16) 

then ( |5.12| ) is still correct for the transfer matrix T{u) = T^^'^\u) of the eight- vertex 
model with diagonal K matrices, which are given by setting fi± = 0. Thus the ansatz 
solves the functional relations of the eight-vertex model. Moreover it is reduces to the 
solution of the six- vertex model as the elliptic nome p ^ 0. These facts, following 
Reshetikhin ||31[, support that the above ansatz could indeed be the eigenvalues of the 
eight-vertex transfer matrix. Therefore the thermodynamics of the models could be 
determined by solving the Bethe ansatz solutions ( [5. 41 ) and ( |5.11| ). 



5.3. RSOS model with open boundaries: In the functional relations for the 
restricted SOS model or ABF model with open boundaries have been constructed. Here 
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we have shown the functional relations of the unrestricted SOS model and the corre- 
spondence between vertex K- and face i^'-matrices. But we have not given the explicit 
form of the face K matrices by this correspondence. This however is a very interesting 
problem. We may obtain the face K matrices of the ABF model by restricting the SOS 
model. Also whether or not the face K matrices following from this correspondence are 



the same as the face K matrices presented in [|T^ is an open question and is presently 
under investigation. It is obvious that the same idea can be applied to higher rank JMO 
IRF models and Belavin's Z„ vertex models. 

The definition of the SOS transfer matrix is not unique. The definition given in 
( [4.15| ) is different from the one presented in because it seems that we cannot get one 
from another one by using simple symmetries like crossing symmetry of the face weights. 
However, the definition given in ( |4.15 ) follows directly from the original formulation given 
by Sklyanin and is good to study a square lattice rotated by 45° (see ll^]), which is 
the natural geometry for an integrable loop model with open boundaries in statistical 
mechanics. This further study will be helpful to generalize the study of Baxter's eight- 
vertex model with periodic boundaries in to the model with open boundaries and to 



well understand the critical behaviour of the face models 33, 34 . 
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Appendix A: Transfer matrices and their graphs 

Here we prove the commuting relation ( |1.22|) for the fused models graphically following 
the Sklyanin's arguments in the study of the unfused six- vertex model 0. For the sake 
of clarity let us represent the R- and i^'-matrices by the graphs. 




(A.l) 
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Therefore we can represent the Yang-Baxter equation ( |1.1| ) graphically by 




V u. 




or 




the reflection equations ( |1.12D and ( |1.13D by 




and 



for (1.12) 

and the unitary and crossing unitary relations by 

1 — 

2 — 





U V 





for ( CT ) 



p{u) 



(A.2) 



(A.3) 



(A.4) 



1- 



-u-2\; 





(A.5) 



Using the graphs the fused R matrix Rip^q) (u) (|1.7|) is given by the elements 





I 


. p 




^ u 






Q 



tp- 



h h 



£2. 



E 



^2- 



ki k2 



1- pi^ 



U22 121 



Ul2 ill 



3p 



-n 



(A.6) 



and the fused K matrices K^\u) ( |1.14[ ) and K_^_ (u) ( |1.15| ) are given, respectively, by 
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the elements 



and 



u 



P{u\q) J2 




u-\-qX—\ 



(A.7) 




E 



u+X 




where Ur^s = u + r\ — s\ and sum over j, 1 is over all possible spins ji, 72, ■ ■ ■ ,jp and 
h,l2,--- Jq satisfying 

p 1 



The spins i and k are 



is and = ^ /c^ 



The indices i, j = —\p, —\p + 1, ■ ■ ■ ,\p and k, I = —^q, + 1, ■ ■ ■ , ^q. 



The fused Yang-Baxter equation and reflection equations (|1.17|) - (|1.18|) are ex- 
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pressed by graphs 





or 



P 




U _ V 




p 



(A.9) 




"^^^ and 



^ (A.IO) 
for (1.17) for ( CT ) 

The unitary (|1.9|) and crossing unitary conditions (|1.1CI| ) for the fused R matrix become 




P' 





(A.ll) 



P' 

q- 





(A.12) 



To prove ([Ol consider T'^P'''\u) T^P'^\ 




(A.13) 



Inserting the operators Y^ and Yj^ respectively into the positions a and b and using the 
crossing unitary condition ( A.12 ) and then using the Yang-Baxter equation (|A.9|) to push 
the vertex with spectra parameter u + v — X + bX through to the most right, we are able 
to obtain 



Pq,b{u+v-X-iii\) 



c 










/ a, 












-lA 








\ V 










u 











(A.14) 
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Then inserting the operators and again into the positions a and c respectively 
and employing the unitary condition (|A.11| ) and then using the Yang-Baxter equation 
( |A.9|) to push the vertex with spectra parameter v — u through again to the most right, 
we have 



Pq,b{u+V-X+b\)pb,q{u-v) 























u-X-lA 








v+a- 


V 




» 






u 











(A.15) 



By the reflection equations ( |A.1CI|) we obtain 



Pq,b{u+V-Mil>^)Pb,qiu-v) Kf{t 

















u-\-qA—\ 




s 


'i-\-lA 


















X: 


u 











K'f\u) 



(A.16) 



Using the Yang-Baxter equation to push the vertex with spectral parameter v — u back 
through to the left and removing the vertices with spectral parameters v — u and u — v 
by the unitary condition ( |A.11|) , then pushing the vertex with spectral parameter v + 
M — A + 6A through back to the left and absorbing the vertices with spectral parameters 
—V — u — X — bX and u + v + bX by the crossing unitary condition ([A.12|) , we flnally obtain 



( |A.13| ) with u and v exchanged. So we have shown that 
required. 



u 



as 



Appendix B: Functional equations and their graphs 

In this appendix we prove the functional equations expressed in Theorem 1. Let us 
consider T^,"^ 



K+{u+q\) 











































u 











K-{u+qX) 



fB.r 



Inserting the operator Y^ into the position a and the identity operator into the position 
b and then using the unitary condition (|A.11|) and the Yang-Baxter equation ( |A.9|) , we 
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are able to obtain 



Pq,i{2u+qX) 



K+{u+q\) 
X 



r 


























u+q\ 






u 









K-(«+<A) 



(B.2) 



It is easy to see that the equation can be split into the sum of the following two parts, 













/ 


U+l 


\ * 




u+q. 


u 




u 







(B.3) 



and 













/ 


u+t 






u+q. 


u 




u 







r2\ 



(B.4) 



where the w-dependent function fq^i{u) is generated from the antisymmetric fusion of 
the model. The underlying models are su(2) type and the antisymmetric fusion gives 
only one independent nonzero element for the R or K matrices. Therefore f^{u) can be 
factorized as 

J^(u) — Oi)^(u + qX — X)oi)+(u + qX — X)(l^(u + qX — X) 

X(l)l(u + qX- X)I+{u + qX- X\q)I-{u + qX - X\q) . (B.5) 

For the six- and eight- vertex models we have shown 

I+{u\q)I_{u\q) = I+{u\l)I_{u\l) = p^1(2m + A) . (B.6) 

Therefore f^{u) depends on q only through the spectra parameter shift u + qX — X. We 
suppress the subscript q and define f^(u) — J^(u — qX + X). Graphically f^(u) is given 

by 



1 



K+{u + X) 



Pi,i(2k+A) k+{u) 



















u+. 










u 














I 











{u + A) 



(B.7) 
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or 



riu) = uj4u)Lu4u)<l>liu)<l>^4u)l4u\q)l4u\q) 



(B.8) 



where 



K+{u+\) 



K+iu) 




k I 
UJ-{u) = 


















u+ 


\ 


1 1 






u 







3 I 



I4u\q) 



j I 





(m+A) 




in) 




-1 




((3-4)A-2u"' 




,((3-3)A-2u 



(B.9) 



(B.IO) 



^ xp-J(2(i-gA+2X) 



(B.ll) 



where the cross with a circle plays the role of antisymmetric operator, 

_ / A; 




The theorem 1 can be seen from ( [B.3| ), ( |B.4| ) and ( |B.8|) 



Appendix C: Intertwining relations and their graphs 

Acoording to the graphs presented in Appendix B and Appendix C the intertwining 
relations between the vertex and face models (f4.3|) can be represented as 





c 








c 




u—v 




i 




V 




j 


b 




17^ 



(C.l) 



where the solid dots mean summation (as do other solid dots in this section). The 
Star-triangle relation (Yang-Baxter equation) of these face weights is graphically 



!■(—* < v-u\c — f( v-u 




(C.2) 
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For the boundaries the correspondence 

'a 



) are given by 





and 




5(c) 
S{h) 




The unitary conditions of the face weights ([4.6|) is 

d d 




a' = p{u)6aa' 



b b 



and the crossing unitary condition ( [4 .71) is 

d d 




, S{a)S{c) 



(C.3) 



(C.4) 



(C.5) 



(C.6) 



The reflection equations for the SOS model ( fl.l4| ) and ([4.13|) can be presented graph- 
ically by 



and 





S{f)S{a) 
S{d)S{g) /' 




S{f)S{a) 
Sib)S{g) f 




(C.7) 



(Ci 
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The transfer matrix of the SOS model with the open boundaries ( [4.15D is presented 
graphically by 



(a|V(M)|6) =K+{i 



bo bo 



U—Vl 


U—V'^ 


— 




■u— llJV 













bN 6iv 



where the boundaries are given by 

Sih) 



■a 

K+{ b 



u :-- 



S{a] 



^^Au + X) K+ (M)^b,e(-M - A) =K+{u)<}b 



and 



Kih 



(C.9) 



(C.IO) 



(C.ll) 



Let us consider V{u) V{u + A) graphically, 



u+X 






U+X 


u+X 






u+X 


u 






u , 


u 






u ■ 



Inserting the crossing unitary condition ( ^.7| ) into above term we have 





(C.12) 



where we have not written down the ratio of the functions S, which come 
crossing unitary condition and we have to take care of them. The sum over 
divided into the symmetric and antisymmetric ones for c = d, 



(C.13) 

from the 
a can be 



1, c ^ c, c + 1, c 



c, c — 1, c c, c + 1, c 



S{c — 1) c,c — l,c c,c + l,c 
\ S{c + 1) c, c + 1, c c,c + l,c 

c,c—l,c S{c — 1) c,c — l,c 
c,c—l,c \ S{c+ 1) c,c+ l,c 
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Therefore ( |(J.13| ) is divided into two terms, which are 



J29sic,b,d) 

b 





where a = c+ l for c = d and a = {c + d)/2 otherwise and 





where c = d and b = c — 1. The factors Qg and Qa are given by 



9s{a,b, c) 



ga{a,b, c) 



S{a-1] 



b = a — l,a = c 



1, otherwise 



S{a-1) 
^5(a + l)' 



b = a + l,a = c 



otherwise 



(C.14) 



(C.15) 



(C.16) 



(C.17) 



The first term is the fused transfer matrix of the model at level 2, which is expressed 
by V^'^\u) = V''q\ The second term is the antisymmetric fusion of the model, which 
gives exactly the factor /q = f^{u). Note that the function f^{u) is dependent the face 
boundary K± matrices. So we finally have the functional equations ( ^4.17|) . 



References 

[1] I. Cherednik, Theor. Mat. Fiz. 61(1984)977. 
[2] E. K. Sklyanin, J. Phys. A21 (1988) 2375. 

[3] H. J. de Vega and A. Gonzalez Ruiz, J. Phys. A26 (1993) L519. 
[4] B. Y. Hou and R. H. Yue, Phys. Lett. A183 (1993) 169. 
[5] T. Inami and H. Konno, J. Phys. A27 (1994) L913. 



26 



[6] H. J. de Vega and A. Gonzalez Ruiz, Nucl. Phys. B417 (1994) 553; Phys. Lett. 
B332 (1994) 123. 



[7; 



[9 
[10 

[11 

[12 
[13 
[14 
[15 
[16 



B. -Y. Hon, K.-J. Shi, H. Fan and Z. X. Yang, Commun. Theor. Phys. 23(1995 163. 

F. C. Alcaraz, M. N. Barber, M. T. Batchelor, R. J. Baxter and G. R. W. Quispel, 
J. Phys. A20 (1987) 6397. 

L. Mezincescu and R. I. Nepomechie, Nucl. Phys. B372 (1992) 579. 

C. M. Yung and M. T. Batchelor, Nucl. Phys. B435 (1995) 430. 
A. Foerster and M. Karowski, Nucl. Phys. B408 (1993) 512. 

A. Gonzalez Ruiz, preprint F.T./U.C.M.-94/1; |hep-th/940lTT8| . 

L. Mezincescu, R. I. Nepomechie and V. Rittenberg, Phys. Lett. B147 (1990) 70. 

S. Artz, L. Mezincescu and R. I. Nepomechie, Bonn-Th-94-19, UMTG-178. 

L. Mezincescu and R. I. Nepomechie, J. Phys. A25 (1992) 2533. 

M. Jimbo, R. Kedem, T. Kojima, H. Konno and T. Miwa, Nucl Phys. B441 (1995) 
437. 



[17] P. A. Pearce and R. E. Behrend, Commuting row transfer matrices and functional 
equations for lattice spin models with hxed boundary conditions. 
After this paper was submitted, two relavant preprints appeared in the HEP-TH 
bulletin- board: 

P. P. Kulish "Yang-Baxter and reflection equations in integrable models" [hep 
th/9507070| ; 



R. E. Behrend, P. A. Pearce aand D. L. O'Brien "Interaction -Round -a - Face 
models with fixed boundary conditions: the ABF fusion hierarchy" [hep-th/9507118. 



[18] P. P. Kulish, N. Yu. Reshetikhin and E. K. Sklyanin, Lett. Math. Phys. 5 (1981) 
393. 

[19] A. N. Kirillov and N. Yu. Reshetikhin, J. Phys. A20 (1987) 1565. 

[20] V. V. Bazhanov and N. Yu. Reshetikhin, Int. J. Mod. Phys. B4 (1989) 115. 

[21] A. Kliimper and P. A. Pearce, Physica A 183 (1992) 304. 

[22] A. Kuniba, T. Nakanishi and J. Suzuki, Int. J. Mod. Phys. A9(1994) 5267. 

[23] Y. K. Zhou and P. A. Pearce, Nucl. Phys. B (1995). 

27 



[24] Y. K. Zhou, Nucl. Phys. B (1995). 

[25] A. Kuniba and T. Nakanishi, Int. J. Mod. Phys. (Proc. Suppl.) A3(1993) 419. 
[26] R. J. Baxter, J. Stat. Phys. 28 (1982) 1. 
[27] R. J. Baxter, Ann. Phys. 76 (1973) 1; 25; 48. 

[28] R. J. Baxter, Phys. Rev. Lett. 26 (1971) 832; Ann. Phys. 70 (1972) 193. 
[29] M. Jimbo, T. Miwa and M. Okado, Lett. Math. Phys. 14 (1987) 123. 
[30] S. Ghoshal and A. Zamolodchikov, Int. J. Mod. Phys. A 21 (1994) 3841. 
[31] N. Yu Reshetikhin, Sov. Phys. JETP. 57 (1983) 691. 

[32] In ||T^ the reflection equations do not follow from the orignal formulation given by 
Sklyanin in [0 and thus they are set up differently. 

[33] J. Cardy, Nucl. Phys. B324 (1989) 581. 

[34] H. Saleur and Bauer, Nucl Phys. B320 (1989) 591. 



28 



